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Effect of electric field on spin splitting in SiGe quantum wells (QWs) has been studied theoretically. 
Microscopical calculations of valley and spin splittings are performed in the efficient sp 3 d 5 s* tight- 
binding model. In accordance with the symmetry considerations, the electric field not only modifies 
the interface induced spin splitting but gives rise to a Rashba-like contribution to the effective two- 
dimensional electron Hamiltonian. Both the valley and spin splittings oscillate as a function of the 
QW width due to inter-valley reflection of the electron wave off the interfaces. The oscillations 
of splitting are suppressed in rather low electric fields. The tight-binding calculations have been 
analyzed by using the generalized envelope function approximation extended to asymmetrical QWs. 
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I. INTRODUCTION 



Understanding the details of semiconductor hetero- 
structure electronic properties has been a key at each 
stage of their applications in the field of information 
and communication technologies. Presently, there is a 
broad interest in spin-dependent properties because they 
have a potential for novel "spintronic" devices, and be- 
yond, because they govern for a large part the possible 
development of semiconductor-based quantum informa- 
tion processing. Spin splitting of electron dispersion rela- 
tions arises from the combination of spin-orbit coupling 
and inversion asymmetry. Besides the contributions of 
bulk inversion asymmetry (BIA) first discussed by Dres- 
sclhaus in 1954 and the structure inversion asymmetry 
(SIA) introduced by Rashba^ the existence of a con- 
tribution due to the breakdown of roto-inversion sym- 
metry at an interface between two semiconductors was 
first suggested by Vervoort et al. in the late 90s<^ 
This "interface inversion asymmetry" (HA) term was 
further documented both by group-theoretical analysis^ 
envelope-function calculations^ and measurements of cir- 
cular polarization relaxation in quantum wells (QWs) 
based on various III-V semiconductors However, in 
these cases, the interface contribution appears in combi- 
nation with BIA and SIA. Pure IIA can exist alone or to- 
gether with SIA in heterostructures of centro-symmetric 
semiconductors like Si-Ge QWs. Previous works^ have 
established the symmetry properties specific to this sys- 
tem where electrons lie in states originating from the bulk 
X z valleys. A general feature of zone-edge conduction 
valleys in bulk materials is their degeneracy (6 for X- 
vallcys, 4 for L-valleys) that gives rise to strong valley- 
coupling when they are folded onto the two-dimensional 
Brillouin zone of a QW^^^^. Valley coupling 
is another manifestation of the local, three-dimensional 
variations of crystal potential at semiconductor inter- 
faces and quantitative estimates require atomistic infor- 
mation which is not available within the fc ■ p theoret- 



ical framework. Parameters describing valley coupling 
must be extracted from microscopic approaches such as 
ab-initio calculations or modclizations using empirical 
paramctrizations like atomistic pseudopotcntials or tight- 
binding approach. The valley coupling strongly depends 
on the crystalline growth direction and shows an oscillat- 
ing behaviour as a function of the number of monolayers 
forming the X z -vallcy or L-valley quantum well. It also 
depends on the overall symmetry of the quantum well and 
for this reason, it can be modified by an external electric 
field. The spin-splittings of in-plane dispersion relations 
in such systems results from the interplay of valley cou- 
pling and spin-dependent terms in the electron Hamilto- 
nian. The case of L-valley QWs formed in the GaSb- AlSb 
system and grown along the [001] direction was first dis- 
cussed in some details by Jancu et alJ^ In that case, 
the leading terms come from BIA invariants specific to 
the L valleys in combination with the L-valley coupling. 
For Si-Ge (OOl)-grown QWs the interplay of IIA with 
X z valley coupling was examined from the fc • p theory 
point of view and semi-quantitative estimates were dis- 
cussed in the frame of tight-binding calculations based 
on the sp 3 s* model.— However, it is well known that 
this simple model cannot reproduce quantitatively the 
properties of zone edge valleys such as effective masses 
and dipolc matrix element. This difficulty was solved in 
the late nineties by the introduction of the extended basis 
sp 3 d 5 s* tight binding modelJ£ More recently, progress in 
the parametrization of this model have led to essentially 
perfect description of the electronic properties of bulk 
Geii and SiJ^ In this work, we use the advanced tight- 
binding model in combination with the envelope func- 
tion approach and calculate the conduction band spin 
splitting resulting from the interplay of valley coupling 
with the IIA and electric-field effects in Si/SiGe quan- 
tum wells. 
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II. POINT-GROUP SYMMETRY ANALYSIS 

In the virtual-crystal approximation for SiGe alloys, 
an ideal (OOl)-grown SiGe/Si/SiGe QW structure with 
an odd number N of Si atomic planes has the point- 
group symmetry T>2d which allows the spin-dependent 
term a{a x k x — <J y k y ) in the electron effective Hamilto- 
nian, where u Xl a y are the Pauli spin matrices, k is the 
two-dimensional wave vector with the in-plane compo- 
nents k x ,k y , and x || [100],?/ || [010]. The QW structures 
with even TV have the T>2h point symmetry containing 
the space-inversion center, the constant a is zero and the 
two-dimensional electronic states are doubly degenerate. 
Under an electric field F = (0,0, F z ) applied along the 
growth direction z the symmetry of QW structures with 
both odd and even numbers of Si monoatomic layers re- 
duces to the point group and the spin-dependent 
linear-fc Hamiltonian becomes 

H^\k) = a(a x k x - a y k y ) + (3(<j x k y - a y k x ) , (1) 

where the second contribution is usually called the 
Rashba term (or the SI A term). In order to establish 
the parity of a, j3 with respect to inversion of the electric 
field we note that, in the T>2d group, the combination 
h(k) — a x k x — cr v ky as well as even powers of F z are in- 
variants while both the combination h'(k) = a x k y — a y k x 
and odd powers of F z transform accor ding to the same 
representation B2 (in notations of Ref. Il9h . Therefore, 
for structures with odd N, the coefficients a and (3 are, 
respectively, even and odd functions of F z . They can be 
presented as 

a(F z ; odd TV) = a + c^F 2 + c^F* + ... , (2) 
/3(F,;oddTV) = cf ) F z +cf ) F 3 + ..., 

where ao = a(0) and <^ n \ch n are field-independent 
coefficients. Similarly, for structures with even TV, the 
linear-in-fc spin-dependent Hamiltonian can be presented 
in the form 

H {1) (F z ;even TV) = F z [dh(k) + C 2 h'{k)} , (3) 

where C\,Ci are even functions of F z . The above rep- 
resentation follows immediately if we take into account 
that, with respect to operations of the D 2 h group, both 
h(k) and h'(k) transform in the same way as the compo- 
nent F z does. 

The aim of this work is to calculate and analyze the 
electric field dependencies of a and j3. For this purpose 
we use the precise nearest-neighbor sp 3 d 5 s* tight-binding 
model 1 ^ and calculate valley and spin splittings in sym- 
metrical QWs in the absence and presence of an external 
electric field. 

III. TIGHT-BINDING MODEL 

To calculate electron subband splittings, we use the 
sp 3 d 5 s* tight-binding theory elaborated by Jancu et alJ£ 



It perfectly reproduces band structure of indirect bulk 
semiconductors as well as electron effective masses, etc. 
In particular the parametrization used in this work re- 
produces the value fco=85% of the conduction band min- 
imum in Si, which was considered as a challenged One 
of the main advantages of this method is a very straight- 
forward treatment of nanostructures. 

In Ref. [Til , we estimated the electron spin splitting in 
symmetrical SiGe QWs using a less detailed tight-binding 
model, namely, the sp 3 s* model, which allowed us to 
understand the main qualitative features of spin splitting 
as well as to demonstrate the possible observability of this 
effect in Si/SiGe heterostructures. 

In the tight-binding model the electron wave function 
is written as a linear combination of atomic orbitals 2 ^ 

\i>,r)=^2C n , v \9 v ,r-r n ), (4) 

where n enumerates atoms in the structure, v runs 
through the set of spinor orbitals at the nth atom. In 
the sp d s* model, this set includes the orbitals s, p v 
(V = x,y,z), d^ (£ = yz,xz,yz,x 2 - y 2 ,2z 2 - x 2 - y 2 ) 
and s* multiplied by the spinors j and J,. We assume 
the basic orbital functions to be orthogonal. Anyway, 
this can be achieved using the Lowdin orthogonalization 
procedure^ Thus, the tight-binding Hamiltonian is pre- 
sented as a multicomponcnt matrix and the Schrodingcr 
equation as an eigenvalue problem 

(Vn,u\ H |tf n ,,„) C n ,y = EC n , v , (5) 

n' m' 

where l^n.i/) = l^vjf — r n ). The Hamiltonian matrix 
elements depend on the relative position of atoms, r n — 
r n i , and chemical type of atoms n and n'. We use here 
the nearest neighbour approximation where the matrix 
elements differ from zero only for neighbouring atoms. 
The detailed procedure of constructing the tight-binding 
Hamiltonian can be found in Ref. |2lJ. Strain effects can 
be included by scaling the matrix elements with respect 
to the bond-angle distortions and bond-length changes^ 

For the SiGe alloy we use the virtual crystal approx- 
imation (VCA) in order to concentrate on the intrinsic 
structure symmetry thus neglecting all effects of disorder. 
Tight-binding parameters were optimized to carefully re- 
produce alloy band structure. We treat strain in two 
independent ways: first, atomic positions used in calcu- 
lations are chosen by using Van de Walle's modeld We 
have also applied Keating's Valence Force Field (VFF) 
models with the SiGe parameters from Ref. [26| and 
found no difference between the continuous and atom- 
istic approaches. 

In addition to the strain dependence of tight-binding 
parameters we have corrected the structure potential 
(see below) with respect to experimentally observed 
Si/Sii-zGc-r conduction-band offseti 27 i 28 

We treat an electric field in the tight-binding approach 
in the following wayd the diagonal matrix elements of 
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the tight-binding Hamiltonian arc shifted due to the po- 
tential of the applied electric field 



(*„,„| H\V n ,, u ,) 

= (tf n>t/ | H\V n ,, v ,) u=Q + U(r n ) 5 nn > <W , (6) 

where U(r) is the electric potential energy. 

Since we are interested in the in-plane dispersion of 
free electrons in a heterostructure we impose periodical 
boundary conditions in the interface plane (001). Be- 
cause of the periodicity in the [100] and [010] directions 
we can introduce the in-plane wave vector k and, for a 
given value of k, construct the tight-binding Hamiltonian 
with a discrete spectrum. For the sake of numerical sim- 
plicity we also use periodic boundary conditions along 
the growth direction [001] taking the barrier layers thick 
enough to exclude the influence of their thickness on the 
calculated values of a and (3. It follows immediately from 
the band structure of silicon and SiGe/Si/SiGe structure 
potential that, neglecting the valley splitting, the elec- 
tronic states with k x = k y = are four-fold degener- 
ate. We will focus on the dispersion of the lowest con- 
duction subband el. The interface-induced valley mix- 
ing leads to a splitting of the state \el,k = 0) into two 
spin-degenerate states denoted + (upper subband) and 
— (lower subband). At nonzero k each subband, + and 
— , undergoes the spin-orbit splitting described by Eq. ([1]) 
with the coefficients a± and (3± for the valley-orbit split 
subbands ±. It is instructive to rewrite Eq. (JTJ) in the 
coordinate frame x' | [ll0],j/' || [110] as follows 



H {1) (k) = (a ± +[3±)<T x ,k y , 



(a± - f3±)<Jy'k x , 



(7) 



Let us introduce the energy difference Ait ] (k || [110]) for 
the states |±,fc || [110]) with the spin polarized parallel 

and antiparallel to [110], and A^fe || [HO]) for the 
states |±,fe j| [110]) with the spin polarized parallel and 

antiparallel to [110]. The modulus of Aio\k) gives the 
spin splitting of the ± subbands and the sign of A^(fe) 
determines the relative position of the split spin sublcvels. 
It follows from Eq. {7J and the definition of Ago (&) that 
the constants a± , [3± can be found from 



,. AL ±] (fc || [110]) + Ag (fc|| [110]) 



,(±) 



B± = lim 



4|fe| 

A^(fe || [110])- A^ } (fc || [110]) 
4|fc| 



Also it should be pointed out that, since the studied 
QWs are quite shallow, the electron dispersion should be 
treated with care. To avoid non-linear effects, very small 
values of k x , k v should be considered. 



IV. RESULTS AND DISCUSSION 

A. Numerical sp 3 d 5 s* model calculations: unbiased 
structure 

In order to test and improve our previous results we 
calculated valley and spin splitting in symmetrical Si 
QW with Sio.75Geo.25 barriers as a function of the QW 
width. For SiGe composition we used optimized tight- 
binding parameters precisely reproducing realistic alloy 
band structure. The strategy for parameterization of the 
Ge-Si alloy in a virtual crystal approximation is as fol- 
low: the parameters of Ge and Si hydrostatically strained 
to the alloy parameter are first calculated and linearly 
interpolated. The small remaining differences with mea- 
sured values are then corrected by fine tuning of a few 
two-center parameters. For interface atoms we use linear 
combination of pure Si parameters and alloy. 

Tight-binding parameters are optimised for bulk ma- 
terials. However, band offsets at the interfaces in the 
heterostructure are also important. For conduction band 
offset we use Shaffler's paper— as a reference. Thus, we 
take a value of 150 meV as a conduction band offset for 
X z valley electrons. 

Figure 1 shows the zero-field results of calculation of 
(a) the valley splitting A v and (b) the constants a± 
for the valley split el subbands in a symmetrical sin- 
gle QW structure with N Si atomic planes sandwiched 
between the thick Sio.75Geo.25 barriers. The splittings as 
a function of N exhibits oscillations, in agreement with 
Refs. mHUll! In Fig. 1 X-shaped crosses depicted 
as vortices of the broken line represent the calculation 
in the envelope-function approximation, sec below. The 
broken line is drawn for the eye. Note that, in order to 
simplify comparison with results obtained by other au- 
thors, Fig. la illustrates the valley splitting not only for 
the QW width region 15^-50 A but also for the region 
60-^70 A. 

Results obtained in the framework of the advanced 
sp 3 d 5 s* tight-binding model show considerable difference 
with our previous estimations The valley splitting is 
significantly smaller, its value decreases by a factor of 
3, whereas the spin splitting increases almost six times. 
This difference is not unexpected since a careful tight- 
binding treatment of Si and its compounds is possible 
in the sp 3 d 5 s* model only. In this regard the goal of 
the previous paper— was to demonstrate that the effect 
of spin splitting in macroscopically symmetrical Si quan- 
tum wells is measurable and to reveal the main qualita- 
tive properties of this splitting. 

The previous theoretical values for valley splitting were 
obtained by both tight-bindin g 10 i 15 and pseudopoten- 
tial methods^ Although the first two papers utilize the 
method of calculation similar to that applied here, a 
straightforward comparison is not possible due to differ- 
ent parametrizations of the Sii-^Ge^ alloy, different alloy 
compositions (x = 0.2 in Ref. [if] and 0.3 in Ref. [l5h and 
conduction band offsets used. However, our results are in 
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FIG. 1: (a) Valley splitting A v as a function of N in a 
Sh^Ge^/Si/Sii-^Ge^ (a; = 0.25) QW structure in the ab- 
sence of an electric field. Solid squares and vortices of the 
broken line represent results of calculations using the tight- 
binding method and envelope function approximation, respec- 
tively, (b) Spin-splitting constants a± versus the number 
N of Si monoatomic layers in the same system (odd N are 
taken in account only). Spin splitting is shown by solid and 
open squares (sp 3 d 5 s* tight-binding calculation, a_ and a+ 
respectively) and x-shaped crosses (envelope function approx- 
imation). 



good agreement with the both estimations. For example, 
for a QW containing 64 Si-atomic layers (32 monomolec- 
ular layers, 9 nm) we obtain for the valley splitting ~ 0.11 
mcV , while Refs. [13 and [l5| present the coinciding val- 
ues of ~ 0.2. Our analysis shows that the valley splitting 
is quite sensitive to the SiGc alloy parameters. By using 
the linear combination of Si and Gc tight-binding param- 
eters for the alloy we could reproduce values of the valley 
splitting obtained by Boykin et alJ^ 

Comparison with Ref. [ID is more straightforward. Fig- 
ure 1 in the cited paper shows dependence of the valley 
splitting on the barrier Ge content for a 16 Si-atomic layer 
QW calculated by the empirical pseudopotential method. 
In particular, for the Sio.75Geo.25/Si/Sio.75Geo. 25 QW the 
valley splitting of about 2.5 meV was obtained^ while 
our estimation is 3.7 meV. This is a good agreement tak- 
ing into account that the two values are obtained in two 
completely different approaches for quite narrow QWs 
where interface effects arc extremely important. 

Figure 2 shows the valley and spin splitting constants 
as a function of the conduction band offset for QWs with 
31, 32, 33, 34 and 64 Si atomic layers. The fifth structure 
is taken in order to provide comparison with Refs. H0IIT5I 
In Fig. 2, in addition to the tight-binding calculations, 
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FIG. 2: Valley (a) and spin (b) splittings for the lowest con- 
duction subband versus the conduction band offset, (a) The 
valley splitting calculated for five QWs with N = 31(B), 
32(A), 33(*), 34(T) and 64(») Si monoatomic planes, (b) 
The spin splitting constants a_ (■,♦) and a+ (□,(>) for 31 
(■,□) and 33(4, 0) Si atomic layers. The fitting by using the 
extended envelope function approach is shown by correspond- 
ing lines. 

we present analytical results on valley and spin splitting 
in the framework of the extended envelope function ap- 
proach. A detailed discussion of the analytical treatment 
is given in subsection II V CI Here we only point out an 
excellent agreement between results for the splittings as 
a function of the QW width and satisfactory description 
of the dependence of these splittings on the band offset. 



B. Numerical calculations in the presence of 
electric field 

Figure 3 demonstrates the variation of spin splitting 
constants with the electric field F z for the el valley-split 
subbands. In accordance with symmetry considerations, 
the calculations show that the spin splitting becomes 
anisotropic in QWs with odd numbers of atomic planes 
and appears in QWs with even numbers of atomic planes. 
The variation of valley splitting is very weak and we do 
not present it here. 

In order to determine the coefficients a± and (3± we 
performed the tight-binding calculation of the spin split- 
tings for the electron wave vectors k || [110] and k || [110] 
and then applied Eq. (|8|) to find the constants a and (3 
directly. The electric field is introduced as a shift of di- 
agonal energies in the tight-binding Hamiltonian. In ac- 
cordance with Eq. ©, we choose electrostatic potential 
to be a linear function of z both inside the QW and in 
the barrier areas near the interfaces. The area of the 
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constant electric field is extended inside barriers enough 
to neglect dependence of the splitting on the choice of 
potential profile. Note that the field values in Fig. 3 are 
small enough to avoid the tunnelling of an electron from 
the QW inside the barrier. 

At zero electric field, f3=0 for arbitrary value of N and, 
similarly, a — for even N. In this case the spin splitting 
A so (fe) is independent of the azimuthal angle of the k 
vector. However, with increasing the field the diversity 
in values of a and (3 for QWs with N = 31,32,33,34 
decreases. 

The further discussion of spin splitting as a function of 
electric field is continued in the next section. It suffices to 
note here that within the chosen range of field values, the 
coefficients a± in QWs with odd N are linear functions 
of F£. In contrast, in QWs with even numbers of atomic 
planes where the spin splitting is absent at zero field, a± 
are proportional to F z . We stress that the field- induced 
change of a± (odd N) becomes comparable to the zero- 
field value of a± in quite weak fields F z ~ 4 • 10 4 V/cm. 

In the previous paper— we developed the extended 
envelope function model in order to demonstrate that 
the spin splitting in macroscopically symmetrical QWs 
is fully defined by interfaces. At non-zero electric field 
two mechanisms are possible, namely, the SIA and IIA 
mechanisms. One of the goals of current research is 
to establish the most important term in realistic QWs. 
To reveal carefully the comparative role of two mecha- 
nisms we present analytical treatment of results shown 
in Figs. 2 and 3 in the framework of the envelope func- 
tion approach. 



Here mi and rrit are the longitudinal and transverse elec- 
tron effective mass in the A valley of the bulk material, zl 
and zr are the coordinates of the left- and right-hand- 
side interfaces, the potential energy U(z) is referred to 
the bottom of the conduction band in Si and given by 

U(z) = V9 b (z) - eF z z 

with V being the conduction-band offset, = 1 in 

the SiGe barrier layers and 8b(z) — inside the Si layer. 
The explicit form of the matrices Vl,r obtained by using 
symmetry considerations is presented in Ref . ITU 

In the zero approximation, i.e., neglecting the valley- 
orbit and spin-orbit coupling, Vl^r = 0, the bispinor is 
given by 



ip(z) = 



ci 

C-2 



X(z) 



where c\ , ci are arbitrary z-independent spinors and the 
function xi z ) satisfies the Schrodingcr equation 



Hox(z) = E X {z) 



(13) 



In the following we take into consideration only the lowest 
size-quantized electronic subband el. 

The next step is an allowance for the interface-induced 
spin-independent mixing between the valleys kg and — fco 
described by the matrices 



V R 



A R 
Ar 



A L 
A T * 



C. Extended envelope function approach 

Here we propose an extended envelope function ap- 
proach 1 ^ to describe valley and spin splittings in the pres- 
ence of an external or built-in electric field. The electron 
wave function is written as 

^(r)=e ik rP[i Pl (z)^ ko (r) + ( p 2 (z)^ ko (r)] , (9) 

where ip±k ( r ) = e ±lfc ° z it±fc (r) is the Bloch function 
at the the extremum points ±fco on the line A in the 
Brillouin zone. The spinor envelopes ifx,(p2 in Eq. @ 
are conveniently presented as a four-component bispinor 



<p{z) 



(10) 



The effective Hamiltonian acting on ip(z) is written as 
a 4x4 matrix consisting of the standard zero- approxim- 
ation Hamiltonian 



w - h2 



1 



kx + k y 
dz mi(z) dz mt(z) 



U{z) 



(11) 



which is independent of valley and spin indices, and an 
interface-induced (^-functional perturbation 



W = V L S{z - z L ) + VrS(z - z R ) . 



(12) 



where A R = Ae" ifc ° a , A L = Ae ifc ° Q , a = z R -z L = Na /4 
is the QW width, ao is the microscopic lattice constant, 
and A is a complex coefficient. 

In terms of the envelope x( z ) the matrix element of 
valley mixing can be written as 

(1|W'|2) = |A| [(xl + Xr) cos(fc a - 4>x) 

+ Kx! -Xfl)sin(fc a-^A)] , (14) 

where |A| and (p\ are the modulus and the phase of A, 
XR,l are the values of the envelope x a t the right and 
left interfaces, x(±a/2), respectively. The valley-orbit 
split states have the energy E± = E ± |(1|W|2)|, where 
Eg is the cigen energy of Eq. (|T3|) . so that the splitting 
is equal to 

A V = 2|(1|H'|2)| 

= 2|A|yxi + x R + 2x 2 L x| cos [2(k a - A )] . (15) 

The bispinors for the upper (+) and lower (— ) states are 
given by 



&=,.(*) = ^ 



Cs 

±e -i0M ( 



(16) 
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where c s is spinor | for the electron spin s = 1/2 and J, 
for the electron spin s = —1/2, </>jvr is the phase of the 
matrix element (|14j) . 

The tight-binding calculations show that the spin split- 
ting is much smaller as compared to the valley splitting 
A v . Therefore, the spin splitting can be considered inde- 
pendently for the upper and lower valley-orbit split states 
(|16[) . The corresponding matrix elements are reduced to 

H' ss ,(k;el,±) = M 1Stls , 

± (Re{M ls , 2s 4cos^ M + Im{M ls>2s '}sin^ M ), (17) 

I 



S h(k) + S' h'(k) 
[X*+P* h{k)+ P '* h'(k)] 



and similar equation for Vl with the coefficients interre- 
lated with A, S, S' ,p,p' due to the mirror- rotation opera- 
tion S4 (odd N) or the inversion i (even N) which trans- 
forms the right-hand-side interface into the left-hand-side 
one. Taking into account the relation between coefficients 
ente ring the matrices Vr and Vl and the notations of 
Ref. [ll| we can write the spin Hamiltonians (fl~7|) in the 
form of Eq. (p} , namely, 

«to(fe;el,±) = a±h(k) + /3±h'(k) . (19) 

For the coefficients a± , f3± describing the spin splitting 
of the valley-orbit split subbands, we obtain 

«± - [Xr- (-l) N XL]S±\ P \H a (ct> p ) , (20) 
(3± = (X 2 R-XI)S'T\P'\H P (^). 

Here 

H a (<f>) = Xr cos (k a a - <j> + <£>)-(-l) N X l cos (k a - - *) 
Hp{4>) = \\ cos (k a -<f> + ®)-XL cos ( k o a ~ § ~ $) , 

$ = arg ( x 2 L e i(koa - M + x^e -1 ** 00 -** 3 ) > 

the parameters S, S' describe the intra- valley contribu- 
tions to the interface-induced electron spin mixing, and 
the parameters p = |p|e^ p , p' = |p'|e*V describe the spin- 
dependent inter-valley mixing. Oscillatory dependence 
of the valley and spin splittings on the QW thickness a 
is caused by interference of electron waves arising from 
intcr-vallcy reflection off the left- and right-hand side in- 
terfaces. 



where the subscript indices 1, 2 enumerate the valleys fco, 
— k and s,s' = ±1/2 are the spin indices, the compo- 
nents Mis, is', Mi 3 ,2s' written as 2x2 matrices Mu,M±2 
are related to similar 2x2 matrices Vr ; h, Vl,u, Vr,12, 
Vt.ia by 

A/ n = Xl v l,u + XrVr,u, M12 = x|Vl,i2 + XrVr,i2 ■ 

The spin-dependent contributions to Vr,h, Vl.ii, Vr^2, 
Vt,a2 are linear combinations of h(k) and h'(k), see 
Rei.M- 



[\+ph(k)+p> fc'(fc)]e- ifc »» 1 , . 

S h(k) + S' h'(k) [ ' 

I 

D. Comparison between tight-binding and 
envelope function approach 

In the absence of an electric field, one has x\ = Xri 
$ = for positive and <3> = ir for negative values of 
cos (k a — (j>\), and Eqs. (fl"5|). ([20]) reduce toii 

A v =4xi|A cos (k a - (j) X )\ , 

a± = (3± = for even N, and 

a± = 2 X 2 L [S ± r]\p\ cos (fc a - <j) p )] , (3 ± = 

for odd JV, where 77 = e 1 * = signjeos (fcoa — 0a)}- The 
curves in Fig. 1 are calculated by using the following best- 
fit set of parameters: |A| = 65 meV-A, 4>\ =0.0137r, \p\ 
=4450 meV-A 2 , p = 0.095?r, S = 650 meV-A 2 . 

Although tight-binding-model values of the coefficients 
in the present work are quite different from the previous 
estimations, the envelope function approach proves its 
adequate description of the valley and spin splittings as 
a function of the QW width. The analytical approach 
with merely five parameters perfectly fits the complex 
microscopical calculation. 

Comparison of new results with experimental data of 
Wilamovsky et ali^ shows the better agreement. With 
the necessary correctional the experimental results give 
-0.34 x 10~ 12 eV-A for the spin splitting constant in 
a 120A-thick QW. The more detailed comparison should 
be done with caution since effects of disorder and built- 
in electric fields can have crucial influence. However, the 
coincidence in an order of magnitude shows that our cal- 
culations agree with the available experimental data. 

According to Fig. 2, the description of dependence of 
spin splitting on the band offset is not so perfect in the 
framework of envelope function approach. In fact, the 
variation in band offset results in a complex (obviously 
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non-linear) behaviour of the parameters in the boundary 
conditions at the interfaces. 

Figure 3 is the main result of this work and we dis- 
cuss it in more details. The agreement between the two 
approaches seen in Fig. 3a shows that the extended enve- 
lope function approach catches the physics of spin split- 
ting induced by the applied electric field. Moreover, the 
fact that this agreement takes place with no addition of 
a SIA term shows that in the system under study the 
IIA contribution is dominating. It should be noted that 
Eq. (f2"0)) for a± contains only parameters which can be 
extracted from Fig. 1, and, indeed, Fig. 3a does not con- 
tain any fitting parameters. Additional fitting parame- 
ters used to describe Fig. 3b are as follows: \p'\ =700 
mcV-A 2 , <j) pl = 7T and S' = 70 mcV-A 2 . 

In the high-field limit xlXr/(Xl + Xr) — > so that 
either x\ < Xr or Xl > Xr, an d Eqs. (frK|, ([20]) transfer 
to 

A v = 2|A| max{xl,x|}, 
a± = 2 sign{Ff +1 }[S ± \p\ cos(0 p - <f> x )] max{x!,x|} , 





/3± = 2 sign{Fj [S't\p'\ cos - <f>x)] max{x|, Xr} ■ 

Since one of the interfaces becomes inaccessible to the 
electron the oscillatory behavior vanishes in strong fields. 

It should be stressed that the parity of the coefficients 
a± and (3± following from the above equations com- 
pletely agrees with the general symmetry considerations, 
see Eqs. © and We also note that a monoatomic 
shift as a whole of the QW position in the structure re- 
sults in an inversion of sign of a± while the values of /3± 
remain unchanged. 

V. CONCLUSION 

The sp 3 d 5 s* tight-binding model has been used to cal- 
culate the electron dispersion in heterostructures grown 
from multivalley semiconductors with the diamond lat- 
tice, particularly, in the Si/SiGe structures. The model 
allows one to estimate quantitatively the valley and 
spin splittings of electron states in the quantum-confined 
ground subband as well as the electric-field dependence of 
the spin splitting. In the employed tight-binding model, 
this splitting is mostly determined by the spin-dependent 
mixing at the interfaces. As a result the coefficients de- 
scribing the Drcsselhaus term in unbiased QWs are oscil- 
lating functions of the odd number N of Si monoatomic 
layers. Under an electric field applied along the growth 
axis a non-zero Rashba term appears in QWs with both 
even and odd Si atomic layers. In small fields, the Dres- 
selhaus term is linear in the structures with even N (D2/1 
point group) and quadratic in structures with odd N 
{D2d point group). Thus, in quite low fields about 10 -4 
eV-cm the spin splitting becomes anisotropic and oscil- 
lations as a function of the QW width are suppressed. 



FIG. 3: Spin-splitting constants a± and /3± for the lowest con- 
duction subband versus the electric field F z calculated for four 
QWs with 31, 32, 33 and 34 Si monoatomic planes. Points are 
calculated in the sp 3 d 5 s* tight-binding model. Lines represent 
fitting by using the extended envelope function approach. 

In addition to numerical calculations, an extended en- 
velope function approach is utilized to interpret the re- 
sults of tight-binding calculations. The inclusion of spin- 
dependent reflection of an electronic wave at the interface 
and interface-induced inter-valley mixing permits one to 
describe quite well the numerical dependencies of the 
valley-orbit and spin-orbit splittings upon the number 
of Si atomic planes and the electric field. 
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